Onset of Fokker-Planck dynamics within a Closed Finite Spin System 
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Relaxation according to Fokker-Planck equations is a standard scenario in classical statistical 
mechanics. It is however not obvious how such an equilibration may emerge within a closed, finite 
quantum system. We present an analytical and numerical analysis of a system comprising sixteen 
spins in which spatial inhomogeneities of the magnetization relax approximately in accord with a 
standard Fokker-Planck equation for a Brownian particle in a parabolic potential. 
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Do closed quantum systems show thermal equilibra- 
tion and if so how does this equilibration proceed? Since 
the Schroedinger equation does not feature any attractive 
fixed point and many systems may ultimately be viewed 
as closed quantum systems, this question has been de- 
bated from the early days of quantum mechanics [H-(^- 
Recently this question has regained considerable atten- 
tion, partially due to experimental progress in the field 
of ultra-cold atomic gases, which may provide the possi- 
bility of observing isolated quantum equilibration exper- 
imentally in a well controlled way [g, [7| ■ On the other 
hand theoretical concepts have been and are being devel- 
oped. Some of these approaches focus on the existence of 
an equilibrium for generic quantum systems and analyze 
properties of pure states concluding that an overwhelm- 
ing majority of them features certain "typical" proper- 
ties which are then identified as equilibrium properties 
[H, HI" 12 1 • The idea that energy eigenstates are likely to 
belong to this typical majority as well is now known as 
the eigenstate thermalization hypothesis (ETH) 0, Il3r - 
[ilj ]. Existence of equilibrium in specific models, such as 
interacting spin or interacting bosonic/fermionic model 
is also a subject of current research, cf., e.g., [17H22| . 
Other approaches focus primarily on the thermalization 
process and aim at mapping the dynamics of crucial 
quantities such as certain expectation values, reduced 
density matrices , etc. onto (quantum) master equations, 
sometimes taking the form of Fokker-Planck equations 
(e.g., for the Caldeira Legett model, etc.). This con- 
cept has been developed in great detail for open quantum 
systems (23[, i.e., the scenario of a comparatively small 
system in contact with a much larger system which is 
usually called a bath. However, concrete examples which 
are not of this "system + bath" structure but can never- 
theless be reliably shown to exhibit features of standard 
phenomenological equilibration seem to be rare. Or as 
stated in [24| "It is an intriguing question whether such 
a Fokker-Planck equation can actually also be found for 
a thermalizing and closed many-body quantum system" . 
This recent publication addresses a system of coupled 
spins in which the z-component of the magnetization is 
reported to follow a Fokker-Planck equation (FPE), close 



to equilibrium and with a non-linear rescaling of the mag- 
netization. Furthermore the FPE-generator is found to 
grow linear in time rather than being constant. The work 
at hand is aimed in a similar direction. However, we 
intend to present a finite spin model in which the (un- 
rescaled) magnetization difference between to equal parts 
of the system relaxes according to a FPE with an time in- 
dependent generator and a diffusion term that is almost 
constant w.r.t. the magnetization difference. 
Our spin model may be described as a finite, anisotropic 
Heisenberg spin-ladder of length N/2 = 8. The Hamilto- 
nian reads 

H = Ho + kV 
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where J = 1 is the coupling strength along the beams 
(labeled by a = L, R) and k the coupling strength along 
the rungs. Specifically we investigate k = 0.2 and k — 
0.15. (This model has been suggested in the context of 
relaxation dynamics in [25j j , it is non-integrable at least in 
the sense of not being accessible by a Bethe ansatz). The 
observable of which we intend to analyze the dynamics 
is the magnetization difference between the beams, x: 
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Since the z-component of the total magnetization of the 
entire system S z is conserved we may restrict our analy- 
sis to a corresponding subspace. Here we choose S z —0, 
which means, loosely speaking, half of the spins point 
up, half down. In the following we are going to compare 
the true quantum dynamics of the magnetization differ- 
ence to the dynamics as generated by a naive, stochastic 
model formulated on the level of a master equation. The 
stochastic model is constructed as follows: Assume there 
is a rate at which mutual spin-flips, i.e., simultaneous, 
contrariwise flips of adjacent spins occur. Let these rates 
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be proportional to i.) the probability with which the 
adjacent spins actually point in opposite directions, ii.) 
the square of the coupling constant between the adjacent 
spins. Since the rates for spin-flips within the beams are 
much larger than for spin- flips across the beams, one may 
assume that the spin-flip dynamics within the chains is 
at an approximate local and instantaneous equilibrium 
during the ongoing much slower spin-exchange dynamics 
between the beams. This implies that the probability for 
an, say, up-pointing spin on the, say right beam, is the 
same for every site on the beam, namely just the total 
number of up-pointing spins on the right beam divided 
by the number of beam sites. Of course the respective 
applies to the left beam and down-spins. Furthermore 
due to the above time scale separation the probabilities 
factorize. Following these principles the rates for transi- 
tions between the magnetization difference subspaces X 
are found to be: 



the projector comprising all states that span a subspace 
with a given magnetization difference X. Let P(e,ae) be 
a projector comprising all energy eigenstates from the in- 
terval [E - AE/2, E + AE/2}. Then we chose our initial 
states as 



f 



Px(0) =z P i 



(0,2)-Px-P(0,2) 



Z = Tr{P {0 ,2)PxP {0 ,2)} 



R 



(X-yX±l) 



f 2XY 
2 T ~N~ ) 



(3) 



The corresponding master equation exhibits already in its 
finite, discrete form strong similarities with a FPE for an 
over-damped particle in a parabolic potential: As long as 
the respective probabilities Px are negligible for the ex- 
treme magnetization differences X = ±JV/4, the first and 
second moments X :— J2x XPx, X 2 := J2x X 2 Px both 
relax mono-exponentially with corresponding relaxation 
rates, Ri = 2 7 k 2 , R 2 = 4(1 -1/N)jk 2 . (The FPE would 
yield corresponding rates with R 2 P = 2R FP ). Regard- 
ing large N and X it is adequate to change to a "mag-/ 3 
netization difference density" z — X/N which becomes 
effectively continuous in the limit of large N . Performing 
a truncated Kramer-Moyal expansion w.r.t. z yields 

d t p(t, z) = -d z {{-d z U{z)p) + l -dl(D{z)p) + 0{dl) (4) 



with U(z) = -yn 2 z 2 and D(z) = 7 /i 2 (I/4 + Az 2 )/N. Up 
to second order in the spatial derivative (i.e., in the large 
length scale limit) this is very similar to the standard 
FPE for a Brownian particle in a quadratic potential 
U(z). The only difference is an additional quadratic 
dependence of the diffusion coefficient D(z) on the po- 
sition. This, however, becomes irrelevant for positions 
sufficiently close to the equilibrium position. 
Above the emergence of standard FPE type dynamics 
from the analogous stochastic model has been estab- 
lished. But do the quantum dynamics as generated by 
the Schroedinger equation (setting h = 1 ) indeed follow 
the stochastic model? This will be analyzed in the follow- 
ing. But before embarking upon a theoretical argument 
let us simply illustrate the existence of this similarity by 
some numerical data. 

The first class of initial states of which we are going to 
display the dynamics is constructed as follows: Let Px be 
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We consider initial magnetization differences X = 0, 1, 2, 
and compute the corresponding evolutions of mean and 
variance, i.e., a(t) — (x(t)} and a 2 (t) = {x 2 (t)) — {x(t)} 2 . 
This is simply done by numerically exact diagonalization. 
Fig. Q] shows the evolutions of the probabilities for the 
various magnetization differences, i.e., Px(t) = (Px{t)) 
for the above initial state with X = 1. The solid lines 
represent the evolutions as obtained from ([3]). In order 
to analyze more initial states and see the relation to the 
FPE dynamics in more detail we plot more data. In Fig. 
[5]the a(i)'s are displayed, together with the correspond- 
ing evolutions as calculated from (The time-scale pa- 
rameter 7 is numerically determined from a projection 
operator approach, see below). Obviously there is good 
agreement. In Fig. [3jthe <x(i)'s are displayed, again to- 
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FIG. 1: Evolution of the probabilities of the magnetization 
differences X for initial state pi(0) (see (JS}) and coupling 
strength k — 0.2. Solid lines correspond to the evolution 
as following from the Schroedinger equation, dashed lines to 
the evolution as following from the (discrete) Fokker-Planck 
equation (J3J) 

gether with the corresponding predictions from ([3]). Ob- 
viously there also is a less pronounced but still reasonable 
similarity. 

Further numerical investigations (the display of which 
we omit here for clarity) show that the agreement of the 
FPE, (or ([3])) with quantum dynamics becomes worse if 
X becomes larger (which is in accord with the findings 
in 24|). The agreement also becomes worse for smaller 
models of the same type, i.e., N — 14, 12, .... This encour- 
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FIG. 2: Evolutions of the expectation values of the magne- 
tization difference for initial states px(0) (see (0)) and cou- 
pling strength k = 0.2. Symbols correspond to the quantum 
evolutions, solid lines to the evolutions as following from the 
Fokker-Planck equation. 



ages the guess that the agreement may become better in 
the limit of X/N < 1 but N -> oo 
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FIG. 3: Evolutions of the variances of the magnetization dif- 
ference for initial states px(0) (see |[5j) and coupling strength 
K — 0.2. Symbols correspond to the quantum evolutions, solid 
lines to the evolutions as following from the Fokker-Planck 
equation. 

Now the question arises whether the agreement between 
the quantum evolution and the FPE dynamics is merely 
incidental. Although we are unable to rigorously spec- 
ify the class of systems and observables to which a naive 
FPE description of the above type applies, we provide 
in the following a consideration which indicates that the 
above specific spin system is most likely not just a sin- 
gular example. This consideration is based on the time- 



convolutionless (TCL) projection operator method [23j. 
We may define a projection superoperator V by: 

Vp = P x ^ L , Px = TT{P x p}, d x = Tr{P x } (6) 

Following the TCL projection formalism yields an equa- 
tion of motion for the Px 's which has the form of a master 
equation: 

Py= £ *%&®Px E R x c yPy (7) 

X^Y X^Y 

(This equation technically applies only for initial states 
that fulfill Vp(Q) — /5(0). However, more recent investi- 
gations imply that its applicability may be considerably 
wider [26]). The rates i?y < ^ L (i) are given as expansions 
in the coupling strength k. Readily assessable are the 
leading (second) order contributions which in many well- 
known cases suffice. They read 

RY°x L2 (t) := / C Y , x (t')dt' (8) 
Jo 

i.e., they are integrals over correlation functions which 
are given by 

Cy,xV) = ^Tr{[V(nP Y }[vl0),Px}} (9) 

where V(t') denotes a time dependence w.r.t. the Dirac 
picture, i.e., unitary dynamics as generated only be 
Hq . A straightforward but somewhat lengthy calculation 
yields 

C Y .x(0) = 8y.x±i Rx -? X±1 (10) 

This result is rigorous and implies that the correlation 
functions, from which the TCL-rates are calculated by 
temporal integration, feature initial values that are di- 
rectly proportional to the rates as obtained from the 
naive stochastic description ([3]). This finding is further- 
more independent of the details of the model, i.e., the 
"left" and "right" part of the system could be any Hamil- 
tonian structure, of course the full system must preserves 
the z-component of the magnetization. For (JTUJ) to hold 
the system neither has to be a spin system, it could also 
be a fermionic (tight-binding) system that preserves the 
overall particle number. In this case the observable would 
be the particle number difference between a left and right 
part of the system. However (fTUf does not necessarily 
render a naive description valid. Such a validity only 
results if i.) the true Schroedinger evolution is well ap- 
proximated by a leading order TCL approach based on 
the projector ([6]) and ii.) all the correlation functions 
Cx±i,x(t') as given in ([S]) decay in a similar fashion and 
on a similar time scale independent of X. Otherwise 



4 



the initial proportionality of the correlation functions to 
the naive rates (1101) does not carry over to a propor- 
tionality of the TCL-rates to the naive rates that can 
be made into an equality by appropriately choosing 7. 
Concerning ii.) one finds numerically that for the above 
specific spin model the correlation functions indeed de- 
cay in similar fashion on a similar time scale of r 3 
for \X\ < 2. From numerical integration of the correla- 
tion functions for those small \X\ we find that in order 
to obtain Rjf±i 2 x (t > T ) ~ Rx±\,x one has to choose 
7 = 0.528. For larger \X\ the dynamics of the correlation 
functions becomes more and more irregular (longer relax- 
ation times, persisting oscillations) accordingly the full 
Schroedinger dynamics are numerically found to be no 
longer in agreement with the FPE. Concerning L), it has 
been reported that the matrix that represents V w.r.t. 
to the basis formed by the eigenstates of Ho must fea- 
ture a certain structure in order for a leading order TCL 
approach to yield valid results (This being true regard- 
less of the existence of so-called time scale separation, 
see below). The features of this structure are somewhat 
complex when stated formally 0, 27 1, however they are 
surely met for matrices the elements of which are chosen 
independently at random according to some distribution 
with mean zero. Although the Hamiltonian at hand does 
not contain any random numbers, V shows features of 
a random matrix, and thus give rise to the applicability 
of a leading order TCL approach. To illustrate this we 
display Figs. 14151 Both are meant to give an impression 
of the structure of the matrix block that corresponds to 
the X = — > X = 1 transition. Fig. |4] shows the indi- 
vidual elements of a 50 x 50 sector of the above block 
taken from the center w.r.t. energy. The elements are 
real, zero on average and show no apparent structure, 
i.e. appear random. To visualize the coarse structure 
we average the absolute squares of the matrix elements 
over small sectors of size 0.12 x 0.12 w.r.t. energy (each of 
these sectors still contains w 1.4-10 4 individual elements). 
Obviously these weights are smoothly distributed on the 
coarse scale with a moderate concentration towards tran- 
sitions between similar energies. In this sense the matrix 
is similar to a matrix the elements of which are drawn in- 
dependently at random according to a distribution that 
only depends smoothly on the distance to the diagonal 
elements. 

Regardless of the structure of the coupling matrix the 
coupling strength is of course crucial for the applicability 
of a leading order TCL projection approach. Generally 
projective approaches are considered "weak interaction 
limit approaches" . And indeed in order to get essentially 
time-independent rates a separation of correlation and re- 
laxation time scales is imperative. This separation grows 
when interactions become weaker. In our specific model 
for coupling strength n = 0.2 those time scales differ be a 
factor of Ri 10, i.e., time scale separation is barely imple- 
mented. This suggests even better results of a projective 
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FIG. 4: Visualization of the fine structure of the matrix block 
corresponding to the X — — > X — 1 transition. Displayed 
are the values of 50 x 50 matrix elements from the center of 
the above block w.r.t. energy. Obviously there is no apparent 
systematic pattern, the entries appear random 
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FIG. 5: Visualization of the coarse structure of the matrix 
block corresponding to the X — — s> X — 1 transition. Dis- 
played are the averages over the absolute squares of the matrix 
elements over sectors of size 0.12 x 0.12 w.r.t energy. Both 
axis correspond to energy. Obviously there is a smooth coarse 
structure with a concentration of weight towards transitions 
between states of similar energy. 



approach at weaker couplings. However, other than in 
the context of infinite systems, for finite systems there 
is also a limit in the direction of weaker couplings below 
which results become worse [27, 28 1. Numerics indicate 
that this limit is already reached at k w 0.15 for the 
model at hand as illustrated by Figs. El [7] Apparently 
the deviations of the true dynamics from the naive de- 
scription already become more pronounced than for the 
case of k — 0.2. Thus for the model at hand the range 
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FIG. 6: Evolutions of the expectation values of the magne- 
tization difference for initial states px(0) (see and cou- 
pling strength k = 0.15. Symbols correspond to the quantum 
evolutions, solid lines to the evolutions as following from the 
Fokker-Planck equation. 
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FIG. 7: Evolutions of the variances of the magnetization dif- 
ference for initial states px{0) (see p|) and coupling strength 
K — 0.15. Symbols correspond to the quantum evolutions, 
solid lines to the evolutions as following from the Fokker- 
Planck equation. 



of possible coupling strengths seems to be quite narrow 
(k roughly between 0.15 and 0.2). However, the lower 
limit of the coupling strengths is expected to approach 
zero very quickly if the system size approaches infinity. 
Thus the possible range of coupling strengths is expected 
to increase substantially for "longer spin ladders" . Inves- 
tigations in that direction are currently under way. 

Concerning the applicability of the above naive 
stochastic description to generic systems it should be 
stated that it is a priori not clear whether or not the cou- 
pling V assumes the above "quasi-random" form when 



represented w.r.t. the energy eigenbasis of the decoupled 
parts Hq. And it is relatively simple to name examples 
for which this quasi-random structure does not emerge, 
like, e.g., any two coupled regular lattices filled with non- 
interacting fermions, etc. Nevertheless the quasi-random 
structure of V is in a mathematical sense typical, fur- 
thermore the emergence of quasi-random matrices from 
non-random Hamiltonians has been found in other mod- 
els, e.g., Ill, 3(|. Thus it appears reasonable to expect 



that naive descriptions will apply to a much wider range 
of models then just the specific spin model discussed in 
the paper at hand. Investigations in that direction are 
currently under way. 

An apparent feature that cannot be in accord with 
any TCL approach of the above type ([6]) (that features a 
unique fixed point) is the fact that the final variances in 
Figs. [3l[7]appear to depend (weakly) on the initial state. 
This feature is absent for the expectation values in Figs. 
El [6] These findings may be viewed as a direct conse- 
quence of the fact that the ETH does not perfectly apply 
to the observable x 2 : For any observable A one expects 
in the long run (ip\A(t)\ip) -> £„Hi|n)|(n|-0(O)}| 2 (|n) 
being energy eigenstates) if sufficiently many incommen- 
surate energies are involved [3l[ . This constant value can, 
however, only be fully independent of the initial state if 
the (n|j4|n)'s are independent of n as the ETH claims 
and as is to good accuracy the case for, e.g. , certain 
Hamiltonians comprising random numbers. |13l Il6l |32 |. 
In the case at hand that would require (n|i|n) = and 
(n\x 2 \n) = 1.1 for all eigenstates n from the relevant 
regime. While the former holds exactly true due to sym- 
metry, the latter is only approximately fulfilled. As a 
consequence the mean values truly approach the same 
equilibrium value regardless of the initial state while this 
is only approximately correct for variances. This resem- 
bles the results in [33|, [34 1 , which also find the possibility 
of some memory effects even in non-integrable systems. 
An issue that has received considerable attention in the 
context of open quantum systems and irreversibility is 
the role of the initial state. Projection operator meth- 
ods, refer routinely to specific initial states (such as fac- 
torizing, thermal, etc.) and it is not obvious within the 
framework of projection operator methods what kind of 
dynamics will result from different initial states [23j . The 
initial states we have analyzed so far are factorizing and 
high entropy. Now we turn towards pure, i.e., zero en- 
tropy initial states. We investigate two types, factorizing 
and correlated, the latter here implies entangled. To eval- 
uate the deviation from the expected relaxation behav- 
ior we calculate 6 — / Q 150 \a>Q(t) — a,FPE(t)\dt/150. Here 
ciQ(t) is the mean magnetization difference as resulting 
from the respective initial quantum state, a,FPE(t) is the 
value as calculated from ©. For the factorizing initial 
states we draw pure states at random (uniformly w.r.t. 
the unitary invariant measure) but from the "5 spins- 
up" subspace on the left and the "3 -up" subspace on 
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the right beam and compute the corresponding product 
states. For the entangled initial states we draw states 
from the full X = 1 subspace at random (which are cer- 
tainly entangled Q) For the product and entangled states 
we draw five initial states for each class and average the 
S's. We find: 5 = 0.0154 for the mixed initial state as 
given in ((5|) with X = 1, 5 = 0.0161 for the pure prod- 
uct states and 5 — 0.0162 for the pure entangled states. 
Obviously all those dynamics follow the FPE evolution 
quite closely. This supports the concept of an equilibra- 
tion that proceeds almost independently of the details of 
the initial quantum state as presented in 26] . 
We investigated the dynamics of the magnetization differ- 
ence between the two beams of a specific 16 spin, ladder- 
type Heisenberg model. This dynamics is found to be to 
some extent in accord with a Fokker-Planck equation for 
an over-damped particle in a quadratic potential, thus 
making the model an example for the emergence of stan- 
dard equilibration within finite closed quantum systems. 
An analysis which indicates that this behavior may be 
generic for a wider range of quantum systems is pre- 
sented. Furthermore the equilibration appears to be in 
reasonable accord with the eigenstate thermalization hy- 
pothesis. Standard equilibration is found for mixed but 
also pure, (to some extent random) factorizing and en- 
tangled initial states which agrees with the concept of 
dynamical typicality. 

We thank R. Tumulka, S. Goldstein, J. Lebowitz and J. 
Ummethum for fruitful discussions on this subject. 
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